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DEVELOPMEHT OF A NONLINEAR VORTEX METHOD 
By 

Osama A. Kandil* 


I. INTRODUCTION 

The present sem-annual report includes the progress of the 
research work conducted under this grcint toward the development of 
reliable nonlinear vortex methods for predicting the steady and un- 
steady aerodynamic loads of highly sweptback wings at large angles of 
attack. It also presents eibstracts of the papers, talks, and theses 
produced through this work. This report covers the period from 
October 1, 1980 to March 31, 1931. 

During this period, research work has been concentrated on 
developing the following two methods: 

1. Modified Nonlinear Discrete Vortex (MNDV) method. 

2. Nonlinear Hybrid Vortex (NHV) method. 

The first method is a modified version of the NDV-method where a 
realistic model of the inviscid leading-edge vortex core is intro- 
duced by the principal investigator. The modified method -alleviates 
the problems previously encountered in predicting satisfactory pres- 
sure distributions on highly sweptback wings. Moreover, it enjoys 
a remarkable success in predicting, for the first time, the latest 
experimental data published by Hummel.^ 

Although the old NDV-method was almost cibandoned in predicting 
the flow details and the distributed aerodynamic characteristics of 
of this problem, the modified NDV-method pinpointed 2 ind cured the 
causes of problems encountered with the old method. The preliminary 
published results of the MNDV-method were well received. In this 
regard, the P.I. has received two letters of recognition, the first from 


* Associate Professor, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508. 

^ Hunmel, D. , "On Vortex Formation Over a Slender Wing at Large 
Angles of Incidence," AGJIRD CP-247, January 1979. 


Professor Hunsael of the Institute of Fluid Mechanics of the Technical 
University in Braunschweig# West Germany and the second from Dr. 

Nangia of the Department of Aeronautical Engineering of the University 
of Bristol, England. Copies of these letters are enclosed in Appen- 
dix A. 

The results obtained from the NHV-method for wings having side- 

edge separations have shown that the calculated spanwise variations 

of the load coefficients are in good agreement with the experimental 
2 

data of Scholz. Compcirisons of the results with those of the NDV- 
method have shown that the NHV-method requires less number of vortex 
panels than the latter method for the same accuracy. 

In the next two sections^, specific points of progress and ab- 
stracts of papers are presented for each method. 

II. MODIFIED NONLINEAR DISCRETE VORTEX (MNDV) METHOD 

1. A paper describing this method was presented at the 12th Congress 
of the International Council of the Aeronautical Sciences, Munich, 
Federal Republic of Germany, October 12-17, 1980. The paper is 
titled, "Numerical Prediction of Vortex Cores from the Leading 
and Trailing Edges of Delta Wings." The paper is attached as 
Appendix B. 

2. A preliminary version of this paper was also presented at the 
28th Meeting of the Aerospace Flutter and Dynamics Council, 
Williamsburg, VA, Oct. 1-3, 1980. 

3. More results of this method for 3/4 turn of the leading-edge 
vortex system were presented at the 17th Annual Meeting of the 
Society of Engineering Sciences, Atlanta, Georgia, Dec. 15-17, 
1980. The paper is titled, "Modeling of Vortex Cores and Feeding 
Sheets of Delta Wings using the Nonlinear Discrete Vortex Tech- 
nique." The published abstract is attached as Appendix C. 


Scholz, V.N., "Kraft and Druckverteilungsmessungen an Tragflachen 
Kleiner Streckung," Forscharb. Ing. Wes., No. 16, pp. 85-92. 
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4. Currently the technique is modified for accurate near-field 

calculations by replacing the concentrated vortex segments with 
vortex panels having linear vorticity distributions. The vorti- 
city functions are expressed in terms of the unknown circulations 
of the original concentrated vortex segments. This is opposite 
to what is being done with the NHV-method. For this purpose, 
the recently developed and highly efficient computer code of 
the velocity field will be used. Applications of this final 
version of the MNDV-method to delta wings with various aspect 
ratios and various angles of attack will be performed to study 
the effects of those two parameters on the formation and inter- 
action of the vortex cores. 

These results will be presented as a peirt of the AIAA paper 
No. 81-1263 to be presented at the AIAA 14th Fluid and Plasma 

Dynamics Conference, Palo Alto, California, June 23-25, 1981. 

# 

The paper is titled, "Recent Improvements in the Prediction of 
the Leading and Trailing Edge Vortex Cores of Delta Wings." The 
other part of the paper deals with the viscous modeling of the 
vortex cores. This work is supported under a separate contract 
with the Naval Air Development Center, Warminster, PA. A copy 
of the paper abstract is attached as Appendix D. 

III. NONLINEAR HYBRID VORTEX (NHV) METHOD 

1. The M.S. thesis of Mr. Li-Chuan Chu is completed in December, 1980. 
The thesis is titled, "A Nonlinear Hybrid Vortex Method for Wings 
Having Side-Edge Separations." The abstract of the thesis is 
attached as Appendix E and a copy of the thesis was delivered to 
the technical monitor. Dr. E. Carson Yates, Jr., who also served 
as a member of the thesis committee of Mr. Chu. 

Currently, Mr. Chu is assisting in the development of the NHV- 
method for leading— and side-edge separations for steady-flow 
problems and for unsteady-flow problems due to small oscillations 
of wings around large mean angles of attack. These problems form 
substantial portions ot his Ph.D. dissertation. 
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The results of the M.S. thesis were presented at the 28th Meeting 
of the Aerospace Flutter and Dynamics Council, NASA-Langley 
Research Center, Hampton, VA, Oct. 3, 1980. 

A highly efficient computer code for the exact calculation of 
the velocity field induced by triangular and quadrilateral vortex 
panels having lineeir vorticity distribution has been developed, 
tested, and completed. Work is underway to develop the necessary 
conditions to eliminate the singularity of the normal component 
of the induced velocity along the panel edges It is a logarithmic 
singularity arising from the geometric discretization of the vortex 
sheet. Mr. Thomas Tureaud, M.S. student, is assisting in this 
problem for his M.S. thesis. He will also be working on converting 
the steady-flow computer code from the CYBER-175 computer to the 
CYBER-203 computer. 

The steady-flow computer code for leading- and side-edge separa- 
tions which employs triangular vortex panels on the free-vortex 
sheets is completed. Currently, it is tested for rectangular, 
delta, and clipped-delta wings. 

The formulation of the unsteady-flow problem using the NHV-method 
is completed. The computer program developed in Item 4 will be 
extended and modified for the unsteady-flow problem. 

Due to the amount of work done this year, due to the complexity 
and importcince of the present problems, and due to the natural 
lack of experience of some of the beginning students, we have 
used the computer relatively heavy in the past 6-12 months. 

Certain measures ^lre now taken to closely monitor the computer 
usage. 
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PROF. DR.-INC DIETRICH HUMMEL 

tM INSTITUT FOR STROMUNCSMECHANIK 
TECHNISCHE UNIVERSITAT DRAUNSCHV/EIC 



3^00 BRAUNSCHWEIG DeC . 10, 

• IINiOOfft WtC A 

TILC^ONl «OAAII AAtAAO UNO AAI t4AS 

HI /Ho/80/204 


Or. Osama A. Kandil 

Associate Professor 

Dept, of Mechn. Engineering and 

Mechanics 

Old Dominion University 

Norfolk, VA 23508 
- U.S.A. - 


Subject: Vortex formation behind delta wings 

Reference: Your paper, presented at the 12th ICAS Congress in Munich, 1980 


Dear Dr. Kandil, 

I thank you very much for having sent me a preprint-copy of your paper presente 
at the ICAS-Congress in Munich. I have read your paper with great interest and 
I congratulate you for your excellent results. 

Unfortunately, I was not able to attend the ICAS-Congress because I was abroad. 

I spent about two months in Korea to deliver a lecture series there. Therefore, 

I missed your lecture in Munich as well as many other interesting papers. But 
I heard that there was a good discussion after your talk. 

This is to let you know that we are presently engaged in similar experiments 
on double-delta wings where we get a twin-vortex system on each side of the wing 
In addition, a wing-canard model is under construction which will enable us in 
near future to measure the vortex system of such a configuration. We will 
publish results as soon as possible and we hope that you will perform the 
corresponding calculations using your program. 

With best wishes for a merry Christmas and a happy new year I am 


sincerely yours. 


(Prof. Dr.-Ing. D. Hummel) 
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• UNIVERSITY OF BRISTOL 

DEPARTMENT OF AERONAUTICAL ENGINEERING 



Proftstor of Aeronautleol En(in«erint: 

L. F. CRAOTREE. Ph.Q.. C.En*.. F.R Ao.S. . 

Talophono 24161 


QUEEN'S BUILDING 

UNIVERSITY walk 

BRISTOL 

BS8 ITR 


Professor O.A. Kandil, 

Mechanical Engineering and Mechanics Department, 

Old Dominion University, 

Norfolk, 

Va. 23508 

U.S.A. 20th November, 19£ 


Dear Professor Kandil, 


Vortex Flow Rep 


As promised, I lave pleasure in sending copies of the following reports. You 
might care to simulate some of the flows with vortex sheet models particularly on 
wings with camber at low angles of attack. 


Nangia, R.K. 


Nangia, R.K. 


Nangia, R.K. 


Nangia, R.K. 


Nangia, R.K. 


Nangia, R.K. 


A Study of Slender, Thin, Conically Cambered Wings with Flow 
Separation. 

University of Bristol, Department of Aeronautical Engineering. 
Report RKN/7701 (October 1977). Revised 1980. 

Slender, Thick, Slarp-Edged Conically Cambered Wings with Flow 
Separation. 

University of Bristol, Department of Aeronautical Engineering. 
Report RKN/7801 (June 1978). 

A Study of Slender Conical Thick Wings and Bodies with Variation 
of Flow Separation Point. 

University of Bristol, Department of Aeronautical Engineering. 
Report RKN/7302. 

Secondary Separation Modelling on Slender Wings. 

University of Bristol, Departm.ent of Aeronautical Engineering. 
Report RKN/7901. 

Separated Flow Past Thin Slender Wings of Arbitrary Planform and 
Spanwise Camber. 

University of Bristol, Department of Aeronautical Engineering. 
Report RK.N/7902. 

Leading Edge Devices, Slats on Swept-Back Slender Wings with Flov 
Separation. 

University of Bristol, Department of Aeronautical Engineering. 
Report RKN/8001 (1930). 


With best wishes. 
Yours sincerely. 
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Nunerlcal Prediction of Vortex Cores of the 
Leading and Trailing Edges of Delta Wings* 
Osama A, Kandil* 

Department of Mechanical Engineering and Mechanics 
Old Dominion University, Norfolk, VA 235C3, USA 


ABSTRACT 

Ihe purpose of the present paper is to pre- 
dict the roll-o? of the vortex sheets emanating 
from the leading- and trai ling-edges of delta wings 
with emphasis on the interaction of vortex cores 
beyond the trailing edge. The motivation behind 
the present work is the recent experimental data 
published by Humel. 

The Nonlinear Discrete-Vorter method (NOV- 
method) IS modified and extended to predict the 
leading- and trailing-vortex cores beyond the 
trailing edge. The present model alleviates the 
problems previously encountered in predicting 
satisfactory pressure distributions. This is 
accomplished by lumping the free-vortex lines 
during the iteration procedure. The leading- 
and trai ling-edge cores and their feeding sheets 
are obtained as parts of the solution. 

The numerical results S ow that the NDV- 
irethod is successful in confirming the formation 
of a trailing-edge core with opposite circulation 
and opposite roll-up to those of the leading-edge 
core. This work is a break through in the high 
angle of attack aerodynamics and moreover, it is 
the first mrerical prediction done on this 
problem. 

NOMENCLATURE 

5R wing aspect ratio 

b wing half span 

b(x) local half span 

Cp static pressure 

LCp net surface pressure 

wing root chord 

LE leading edge 

LCC leading-edge core 

TE trailing edge 

TEC trai ling-edge core 

xyt wing-fixed coordinates, origin at wing 

nose, x-axis is along the root chord, 
y-axis is perpendicular to wing 
xyi wind-fixed coordinates, origin at the 

trailing edge, x-axis is parallel to the 
free stream directirr 

v-,v- V- components of v in the wind coordinate 
* * system 

vg- component of v in planes x • constant. 



angle of attack 

diTen$1on1ess_wind-f1xed coordinates, 
i • x/b, . • y/b, ; • i/b 


I DESCRlPnON OE THE FLOW FIELD 

The flow field around highly swept wings at 
moderate to large angles of attack is character- 
ued by flow separation's from the leading and side 

•this research work is sponsored by NASA Langley 
Research Center under NSG lEoO, Or. E.*^Carson 
fates, Jr is the technical monitor 
‘Associate Professor, AIAA member 


edges due to strong cross flows. The flows from 
the upper and lower surfaces of the wing leave at 
these edges forming free-shear layers. The sepa- 
rated free-shear layers roll up spirally and form 
tow vortex cores which are continuously fed by 
the vorticity shed from the attached boundary 
layers on the wing surfaces through the free-shear 
layers. This flow separation is known as the 
"primary separation." It has a dominant effect on 
the aerodynamic charateristics due to the large 
strength of its vortex core and it proximity to 
the upper surface of the wing. It generates a 
large suction pressure peak on the upper surface 

under the primary vortex core.^*^ 

The adverse pressure gradient outboard of the 
suction peaks affects the boundary layer flow on 
the upper surface and “secondary separation" from 
the wing surface occurs. The secondary separated 
flow forms either an additional free-shear layer 
or a bubble depending on the angle of attack and 
the wing aspect ratio In tne range of moderate 
to large angles of attack, the secondary free-shear 
layer rolls up spirally in an opposite sense to 
that of the pn-ary free-shear layer and forms a 
secondary vortex core with a strength much smaller 
than that of the primary core and of opposite 

strength.^*® 

The effect of secondary separation on the 
upper-surface pressure distribution depends on 
the type of boundary-layer flow on the upper 
Surface. For a laminar boundary layer, the secon- 
dary core produces another lower pressure peak 
between the secondary and primary lines of 
separation. For a turbulent boundary layer, such 
a pressure peak is hardly noticed, and the pres- 
sure peak corresponding to the primary separation 
is higher than that of the laminar boundary- layer 

flow.^ Therefore, when inviscid models are used 
,0 predict the pressure distribution, the calcu- 
lated results should be compared with the experi- 
mental data corresponding to a turbulent boundary- 

layer flow ® 

A third type of flow involving a “tertiary" 
separation may occur between the lines of secondary 
and primary sep-ration due to an adverse pressure 
gradient created by the secondary vortex core. 

The free-shear layer emanating from the 
trailing edge of a wing with highly swept-back 
leading edge is of opposite strength to that of 
the primary free-shear layer. Beyond the trailing 
edge, it rolls up spirally In an opposite sense to 
that of the primary free-shear layer and forms a 
vortex core The trailing core has the same sense 
of strength and core rotation as the secondary 
vortex core, although they originate from different 
phenomena. However, their interaction beyond the 

trailing edge is still unknown. Figure I 
(reproduced from Reference 5) shows a schematic of 
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the vortex fomatlon behind e slender delta wing. 
Figure l-b shows that the secondary vortex sheet 
and the tralllng*edge vortex sheet have the same 
sense of rotation and the same spanwlse location 
Their sense of rotation Is opposite to t*iat of 
the primary vortex sheet. 

The purpose of the present paper Is to pre- 
dict the roll up of the free-vortex sheets 
emanating from the leading and trailing edges of 
delta wings. The main emphasis Is to confirm the 

measurements and conclusions of Humsl^ that th« 
traillng-edge sheet rolls uP Into a vortex core 
within a distance of 1/4 root chord from the 
trailing edge, that Its rotation 1$ opposite to 
that of the primary vortex sheet, and that the 
trailing edge core Is of origin different from the 
secondary-vortex core although they have the same 
sense of rotation. 

II. EXISTING MATHSHATICAL HQOEIS 

The description given above shows that the 
flow field contains four different vortex cores, 
primary, secondary, tertiary, and traillng-edge 
vortex cores. The secondary and tertiary cores 
are due to viscous phenomena and cannot be modeled 
by using Invl scid models. However, their effects 
are small partlculary when the turbulent boundary- 
layer flow exists on the upper surface of the wing. 
Hence, they are neglected. 

, In most analyses, the attached boundary 
layers on the wing surfaces are represented by 
bound-vortex sheets while the separated free-shear 
layers are modeled by free-vortex sheets. The 
free-vortex sheets Join the bound-vortex sheets 
along the separation lines which are known a priori 
for wings with sharp edges. Moreover, we assume 
that vortex-breakdowr points are far downstream so 
that the primary core site and Its variation In 
the vicinity of the wing are neglected. In fact, 
this assumption limits the large angles of attack 
at which the Invlscld model Is applicable. 
Furthemore, the flow outside the bound-vortex 
sneet(s) and the free-vortex sheet(s) can be 
assumed Irrotatlcnal. 

Within these assumptions, the resulting 
potential flow model represents the main features 
of the real flow to a high degree of accuracy. 

The literature contains several steady and 
unsteady invlscld-flow models with various degrees 
of limitations and drawbacks. 

The first group of models uses slender-body 

theory and conical flow assumption. These 
models satisfactorily predict the pressure 
distribution over the front portion of the wing 
surface. In the rear portion, the models fal', 
to predict a satisfactory pressure distribution 
because Xutta condition Is violated at the trall- 
'ng edge. Such models are limited to simple delta 
planforms. 

The second group of Models uses a nonlinear 
19 11 

discrete-vortex method. In these podels, the 

bound-vortex sheet and the free-vortex st;eets are 
approximated by a set of concentrated /ortex lines. 
The bound-vortex sheet Is replaced by a bound- 
vortex lattice, while the free-vortex si ret Is 


replaced by segmented free-vortex lines (In the 
case of steady flow) or by a growing free-vortex 
lattice {In the case of unsteady flow). The 
boundary conditions are satisfied at certain con- 
trol points on the bound- and free-vortex system 
using an Iterative technique. Excellent agreemenc 

was found^*’ between calculated and experimental 
total aerodynamic cnaracterl$f*cs, and the agree- 
ment between calculated and experimental section 
characteristics was satisfactory for wings with 
only side-edge separation. For wings with leading- 
edge separations, hc-<ever, the agreement was less 
than satisfactory for sene cases. 

Although the discrete-vortex model has for 
many years worked very well for attached-flow 
32 

problems , when vortex-type separation from 
leading edges and/or tips occurs, the free-vortex 
system lies close to the llftlng-surlace bound 
vortices and results are found to be sensitive 
to variations In the shapes of the quadrilateral 
vortex elements and the relative lengths or the 

vortex segments. The be'-t agreement was ob- 
tained by using the vortex arrangement developed 
In reference 21, 

The third group of models employs doublet- 

panels, In this nethod, the wing and Its 
free-vortex sheets are divided Into networks of 
quadrilateral panels. Each panel of the networks 
representing the wing has a biquadratic local 
doublet distribution and a bilinear local source 
distribution. The panels of networks representing 
the free-vortex sheets have biquadratic local 
doublet distributions. Source and doublet splines 
are used to express the distributions of singu- 
larities on the networks in terms of dlsemte 
values of singularity strength at certain standard 
points on each network. The boundary conditions 
and continuity of singularity strengths across 
abutting networks are enforced at certain standard 
points on each network. The results of this method 

are generally good when the solution converges. 
Apparently, the difficulty In obtaining convergence 
Is due to the failure in satisfying the continuity 
of the derivatives of the doublet strength across 
abutting networks. This Is equivalent to the 
existence of concentrated vortex lines between 
abutting networks. 

The doublet panel method was extended'*^ to 
include the effect of entrainment of the primary 
vortex cores through an empirical approach. The 
results Indicated that the entrainment Increased 
the normal-force coefficient substantially over 
the experimental values. 

An excellent review of the Invlscld models 
discussed above is given in reference 4i. 

Reference 42 evaluates the AGARD S/nposltn on High 
Angle of Attack Aerodynamics, where many relafeJ 
works are discussed. 

The fourth group of moccls employs a nonlinear 

hybrid vortex rethod^^*^^. In this rethon, con- 
tinuous 'ortlclty and ,ortex-llne representations 
of the Ing and Its separated free-shear layers are oScd 
Continuous vortlclty Is useH In the near-field 
calculations while discrete vortex-lines are used 
In the far-fleld calculations. 
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The wing end Its free-shear layers are divided 
into quadrilateral vortex panels having first- 
order vortlcity distributions. The aerodynamic 
boundary conditions and continuity of the vorticity 
distributions are satisfied at certain nodal points 
on the vortex panels. An iterative technique is 
used to satisfy these conditions In order to 
obtain the vorticity distribution and the shapes 
of the free-vortex sheets. The calculated net 
Surface pressure and section nomal-force 
coefficients for wings with tip separation are 
in good agreement with the experimental data and 
the nonlinear discrete-vortex method. Work is 
underway to include wings with leading-edge 
separation and to improve flow modeling by using 
higher-order vorticity distribution and nonplanar 
panels. 

So far, none of the four groups has been used 
to predict the leading- and trailing-edge vortex 
cores. The present work is the first one in which 
such a prediction is considered and successfully 
accomplish this goal . 

in. TH E EXISTING WMINEAR DlSCRETE-tORTEX 
METHOD~ThDV-M£lHOO) 'Af.P lYS bSAwBACkS 

After the AGARD Symposium on High Angle of 

Attack Aerodynamics^, the first attempt to predict 
the vortex coretwas done by the NOV-method which 
was developed in reference 21 and was later re- 
fined in reference 31. Next, a critical evalua- 
tion of the existing model is given in order to 
pinpoint its drawbacks. 

In this mode', as well as in all the other 

19-31 

existing models of the NDV-method no attempt 

has been made to model the leading-edge vortex 
core. The vortex sheets emanating from the lead- 
ing and trailing edges of a delta wing are replaced 
by a system of vortex lines. Each vortex line is 
divided into straight, short vortex segments and 
a semi-infinite vortex line. The dimensionless 
length of each vortex segment is equal to a unit 
length (the chord length of a quad i lateral bound- 
vortex element is taken as the chai actenstic 
length of the model ). 

14 IS 

Starting with an initial guess ’ for the 
inclination of the free-vortex lines representing 
the free-vortex sheets, the boundary condition on 
the bound-vortex lattice representing the bound- 
vortex sheet are satisfied at certain control 
points. The resulting set of linear algebiaic 
equations is solved and the circulation distri- 
bution is obtained. Next, with the knownci rcula- 
tion distribution, the positions of the vortex 
segments of the fr»e-vortex lines are adjusted in 
order to satisfy the boundary conditions on the 
free-vortex sheets. These two steps of calcula- 
tions represent one iterative cycle. Several 
iterative cycles are perfamed until the positions 
of the vortex segnents or the circulation distri- 
bution converge. 

Figure 2 shows a typical convergei} solution 
of the system of free-vortex lines in three views 
for a delta wing of aspect ratio of unity and 15’ 
angle of attack. The^ilan view also shows the 
arrangement of bound-vortex lattice In t^e three 
dimensional view the leading-edge core (LEC) is 
shown. This core is calculated after fe solution 
converges and it represents the centroid of the 


leading-edge vortex system . It is calculated in 
cross-flow planes. 

The comparison of height and spanwise positio 
of the calculated centroid with those of the 

1 4 

measured leading-edge core ' is given in figure 3 
Although this comparison was encouraging, the 
calculated centroid does not model the physical 
vortex core where the vortex jore is continuously 
fed with vorticity from the leading edge through 
the free-shear layers Using the system of free- 
vortex lines, the total -aerodynamic loads were 
accurately calculated but the calculated surface 
pressure distribution was unsatisfactory. • On the 
other hand, using the centroid of the leading-edge 
vortex system instead of the leading-edge vortex 

system, the calculated surface pressure improved^* 
(this does not satisfy the no-penetration condition 
on the wing since thecentroid is calculated after 
the solution converges.) It should be noted here 
that replacing the leading-edge vortex system by 
its centroid is similar to the model used by 

legendre*®. 

Figure 4 shows a recent converged solution of 
the system of free-vortex lines with a long de- 
formed wake. It can be seen that the trace of the 
trailing-edqe vortex system in cross planes indi- 
cates that the sheet tends to deform upwards 
showing a tendency to form a trailing-edge vortex 
core However, the cross-flow planes taken further 
downstream shows that the free-vortex lines leap- 
frog. This does not represent the real flow. 

Figure 5 shows four of this cross-flow planes 
taken perpendicular to the wind direction. 

It IS clearly seen from the few examples 
given above, that the existing model of the NOV- 
nethod does not real istical ly model the leading- 
and trailing-edge vortex cores. Therefore, the 
model and the numerical technique must be modified 
in order to obtain realistic vortex cores. 

IV. MODIFICATION OF THE f.DV-METHOO AND MODELING 
OF THrL'EASrN5:Eofi'£ 

The most obvious drawback of the model is the 
lack of a realistic model of the leading-edge core 
and its feeding vortex sheet. Modeling of the 
leading-edge vortex core and its feeding sheet 

had been first introduced by Brown and Michael 
However, the leading-edge sheet in their model was 
taken as a planar surface and hence it does not 

represent the real flow. Mangier and Smith^^ and 

Smith^^ introduced the first realistic model of 
the vortex core and Its feeding sheet. Their 
model follows the slender body theory and hence it 
does not account for the deformation of the 
leading-edge vortex sheet in the chordwise direc- 
tion The numerical models given in references 
34-37 and 40 also do not account for the deforma- 
tion of the leading-edge vortex sheet in the chord- 
wise direction. Consequently, the leading-edge 
vortex sheet does not feed the leading-edge vortex 
core beyond the wing trailing edge. Moreover, the 
model given in reference 37 represents the free- 
vortex sheets beyond the trailing-edge by using 
a prescribed shape which was called a “f>xed design 
wake" Therefore, none of these models are capable 
of predicting the vortex sheet deformation beyond 
the trailing edge. 
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!n the present modified model of the NOV* 
method, alt these drawbacks are avoided. In tne 
present method, the old model Is used In the 
first two Iterative cycles, this Is the coarse 
model. The centroid, of the teadlng-edge vortex 
system of the coarse model Is calculated and is 
used as the preliminary leading-edge core. Next, 
the leading and tralling-cdge free-vortex lines 
are replaced by finer vortex segments, this Is 
the fine model. Then, the Iterative cycles 
proceed. During each Iterative cycle, each free- 
vortex line Is allowed to rotate three-dimension- 
al ly around the most recently calculated centroid 
for a prescribed number of turns (1/4, 1/2, 3/4, 
or 1 turn). Once this Is accomplished, the re- 
maining vortex segments of this line Is dum.ped 
Into the centroid. This technique Is followed 
until the circulation distribution converges. 

Figure 6 shows typical solutions at different 
Iterative cycles. It can be seen that after one 
Iterative cycle, the system of free-vortex lines 
of the coarse model shows good deformations. The 
converged solution, indicated by ITER • 6, shows the 
leading-edge core and its feeding free-vortex 
lines. It can also be seen that the free-vortex 
lines continues to feed the LEC beyond the , 
trailing edge. The trailing-edge core is also 
indicated on the figure. This will be clearly 
seen In the cross-flow planes discussed In the 
next section. 

Figures 1 and 8 show converged solutions 
using a 12 X 12 bound-vortex lattice and different 
turns of the free-vortex lines. 

V. NlflERICAL EXWLES ANO CCHPARISONS WITH 
HlhWEL"‘"s TxP'fftTMTNTAi:' Wn',"'PRtDl£TION OF THE LEC 
and' TeC 

The calculated circulation of the free-vortex 
lines eininatlng from the wing tralling-edge was 
found to be of opposite sign to that of the free- 
vortex lines emanating from the wing leading-edge. 
This is consistent with Hummel ‘s measurements. 

This Is not a surprise and it can be explained as 
follows: 

The leading-edge vortex core creates large 
suction pressure peak and hence the pressure 
continusouly rises In the sparrwise direction 
from the wing axis to the location of the suction 
peak (this Is completely opposite to the spanwise 
variation of pressure for wings with large aspect 
ratio and small angles of attack; the classical 
llnearlaed wing theory.) Consequently , the 
Circulation of the spanwise bound-vortex segments 
Increases in the spanwise direction. In order to 
satisfy the spatial conservation of circulation at 
the nodes of the bound-vortex lattice, the span- 
wise Increase of circulation requires the circula- 
tion of the chordwise bound-vortex segments to be 
of opposite sign. In the present model, the 
leading-edge vortex system originates from span- 
wise-vortex segnents at the leading edge while 
the trailing-edge vortex system originates from 
chordwise-vortex segments. 

According to the difference in sign of the 
circulation of the trailing-edge vortex system 
from that of the leading-edge vortex system, one 
expects that the trailing-edge vortex system rolls 
up In an opposite sense to the roll-up of the 


leading-edge vortex system. 

Figures 9-13 show comparisons between the 

measured^ and calculated leading- and tralllng- 
eoge sheets and flow directions in cross-flow 
planes perpendicular to the wind direction. The 
calculated results are drawn with the same scale 
as that of Hunnel's measurements. The calculated 
results show a remarkable success of the new 
modal. The predicted slaes and locations of the 
leading- and trailing-edge vortex sheets are 
very close to those of measurements. One can also 
see that the heights and spanwise positions of 
the LEC are successfully predicted. The sequence 
of the calculated cross flow planes shows that 
the roll-up of the trailing-edge sheet occurs In 
an opposite sense to that tta lesdbf-cd^ sheet. 

Figure 13 shows the formation of the trailing- 
edge core (TEC). The difference between the 
pred.cted and measured positions of the TEC is 
due to the number of turns taken In calculating 
the leading-edge vortex system. It Is expected 
to get the best comparison when one full turn is 
considered. With one full turn of the leading- 
edge vortex system, the roll up will tighten and 
larger induced velocities are created which In 
turn win deform the trailing-edge vortex system 
upwards and leftwards as well. This refinement 
Is currently tested. 

Figure 14 shows a converged solution and the 
cross flow planes of a delta wing of aspect ratio 
of 1.45. The formation of the tralllng-edge 
core develops earlier than that of the aspect 
ratio of unity. Work Is underway to test different 
aspect ratios and different angles of attack. 

Figure 15 shows comparisons between the 
predicted and measured static-pressure contours 
In different cross flow planes. The predicted 
sizes, locations, and levels of the pressure con- 
tours are in good agreement with those of the 
measured data. 

Figure 16 shows comparisons of the predicted 
and measured spanwise net surface pressure varia- 
tion at different chordwise stations. With one 
full turn of the leading-edge vortex system, the 
predicted suction pressure peak will reach the 
measured one particuarly In the front portion of 
the wing due to the proximity of the feeding sheet 
to the wing surface. Moreover with one full turn, 
the roll-up tightens and the predicted pressure, 
away from the location of the suction peak, drops 
to match the experimental data. 

The newly developed computer code Is computa- 
tionally efficient. On a CYBER-175, a typical 
solution of the 12 x 12 lattice including the cal- 
culations of five cross-flow planes takes about 
10 minutes of CPU ti-iie. 

CONCLl'OlNg REHARkS 

The newly developed NOV-method Is a break 
through In the area of the high angle of attack 
aerodynan.es It enjoys a remarkable success in 
predicting, for the first time, the latest experi- 
mental data published by Humnel. The present tech- 
nique cures the problems previously encountered with 
this method. The computer code which employs the 
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present technique is highly efficient when com- 
pared with the present computer codes. It is the 
only existing code which is capable of predicting 
the LEC and TEC beyond the trailing edge. Work 
is underway to refine the technique by using one 
full turn of the leading-edge vortex system. 
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Figure S. Leading- and trail tng-edge sheets 
behind a delta wing in planes 
perpendicular to the wind direction, 
AR • 1, a • 17'’, 7x7 lattice. 
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Hgure 16. Spanwise net surface pressure varia- 
tion at different chotdwise stations, AR • 1, 

. ■ 20 S’, 12 » 12 lattice. 1/2 turn. 
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Prediction of the Leading and Trailing Edge 
Vortex C ores of Delta Wings . 0. A. Kandilf 
Associate Professor, Department of Mechanical 
Engineering and Mechanics, Old Dominion Univ., 
Norfolk, VA. 

The recently Modified Nonlinear Discrete Vortex 
Method (MNDV-Mathod) has shown a remarkable 
success in predicting, for the first time, the 
latest experimental data published by Hummel 
on vortex formation over a slender delta wing 
at an angle of attack of 20.5°. 

- This paper presents the recent developments in 
the HfIDV-Method used to accurately predict the 
location of the trail ing-edge vortex core and the 
surface pressure distribution. Also, it presents 
more numerical results of this technique for delta 
wings with various aspect ratios and various 
angles of attack in order to study the effects of 
those two parameters on the formation and inter- 
action of the vortex cores. 

Moreover, a viscous core model, based on the 
boundary-layer-like approximations, is presented. 
The viscous equations, with outer-edge boundary 

co nditions obtained from the inviscid m.odel , are 

““integrated using a finite-difference marching! 
^technique. ' 
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ABSTRACT 



» 

A NONLINEAR HYBRID VORTEX METHOD FOR 
WINGS HAVING SIDE-EDGE SEPARATIONS 

Li-Chuan Chu 

Old Dominion University, 1980 
Chairman; Dr. Osama A. Kandil 

A Nonlinear Hybrid-Vortex method (NHV-method) has been developed 
for predicting the aerodynamic characteristics of wings exhibiting 
side-edge separations. The present method is a coupline between vortex- 
panel and vortex-line methods. In the near-field calculations, vortex 
panels are used while in the far-field calculations, vortex lines are 
used. 

The wing and its free-shear layers are divided into planar 
quadrilateral panels having first-order vorticity distribution. The 
aerodynamic boundary conditions and continuity of the vorticity distri- 
butions are imposed at certain nodal points on the panels. An iterative 
technique is used to satisfy these conditions to obtain the vorticity 
distributions and the wake shape as well. To expedite the convergence 
of the iterative technique and to avoid excessive distortions of the 
free-vortex panels, good initial guesses of the shapes of the free- 
vortex sheets, obtained from the Nonlinear Discrete-Vortex method 
(NDV method), were used. 
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This method is used to calculate the distributed and total steady 
loads on thin flat rectangular wings of different aspect ratios and at 
different angles of attack. The agreement between the calculated 
results and the available experimental data is satisfactory. 


ACKNOWLEDGEMENTS 


It gives me great pleasure to express my profound gratitude to 
Dr. Osama A. Kandil, my advisor and chairman of my committee, for 
constantly providing me his guidance, advice and encouragement during 
the development of this research work and during my education as well. 

I would like to thank Dr. E. Carson Yates, Jr. of the Aero- 
elasticity Branch of the NASA-Langley Research Center and a member of 
my committee for his many discussions and valuable advice during the 
development of this work. 

Thanks are also due to Dr. Gene L. Goglia, Chairman of the 
Mechanical Engineering and Mechanics Department and a member of my 
committee, for his support; Mr. Robert N. Desmarais and Dr. Robert M. 
Bennett, both of the Aeroelasticity Branch of the NASA-Langley Research 
Center, for their valuable help during the development of the computer 
program. 

This work has been supported by the NASA-Langley Research Center 
under grant No. NS6 1560; Dr. E. Carson Yates, Jr. is the technical 
moni tor. 

Finally the wonderful typing job of Ms. Beverly Dorton of the 
Electrical Engineering Department, is highly appreciated. 


31 



f 


End of Document 


/ 



